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1. Introduction 


In recent investigations [2], [3], [9], Banks, Inman, and their colleagues carried out a series 
of experimental and computational investigations on linear damping models for composite 
material flexible structures. Other investigations [1], [7] involving experiments with flexible 
structures point to the need to understand nonlinear aspects (in particular, damping) of 
flexible structures. In this paper, we present a first step toward development of a rigorous 
theoretical foundation for computational methods to study nonlinear mechanisms in such 
structures. 

Traditionally, nonlinearities in structural models usually arise when one needs to con- 
sider large amplitude displacements. However a rather substantial engineering literature 
(see Chap. 7 of [19]) exists on a wide variety of “nonlinear” effects in structures. There 
are numerous sources proposed for these nonlinearities; they can be geometric (nonlinear 
stretching; large curvature) or material (nonlinear stress-strain laws; nonlinear damping) 
in nature. In mathematical models these can be manifested in nonlinear boundary con- 
ditions, nonlinear stiffness operators, and/or nonlinear damping operators. For example, 
in the usual derivations [11], [20], [24] for the Euler-Bernoulli theory of beams or the 
Love-Kirchoff theory of plates, one must make several “linearizing” assumptions. The 
commonly postulated Hooke’s law is a linear stress-strain constitutive law. Even under 
this assumption, one has a nonlinear moment- curvature relationship (and hence a non- 
linear stiffness operator in the resulting partial differential equation model) unless one 
makes the usual approximations (e.g., see Chap. 11 of [20]) in the curvature expression 
1 /R = UxT j [l + a*] 3/2 arising in elementary beam theories. 

While the ideas presented below can be used to develop a methodology to study nonlinear 
stiffness operators (which, of course, can be of great importance in motions involving large 
geometries or nonlinear material properties), we shall focus here on a theory for nonlinear 
damping operators. Our approach will allow one to study damping mechanisms that do not 
decouple across the modes of the structure and this has substantial practical significance. 
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The efforts of Balakrishnan, Taylor, and their colleagues (see [1] and the discussions and 
references therein) have produced evidence that in complex structures such as the SCOLE 
structure at NASA Langley, one cannot use linear damping models where the damping 
decouples across the “linear undamped modes” to successfully describe experimental ob- 
servations. The presentation in [1] suggests that modal-based nonlinear damping might 
be adequate. However, in light of the results in [3], [9], where it is found that for even 
simple structures one can be in the unfortunate situation where damping does not decou- 
ple across the modes, this issue is far from settled. Indeed, it is not clear whether such 
damping should be modeled with a nonlinearity or whether it might be adequately mod- 
eled with a linear mechanism that does not decouple (i.e. cannot be described via modal 
damping ratios in the usual engineering manner). Whatever the case, it is clear that fur- 
ther investigations are necessary. Thus one requires methods to study linear and nonlinear 
damping mechanisms in models of second order systems where one does not assume that 
the damping model can be simplified by modal representations. 

Our goal then is to develop such a methodology for nonlinear systems. The theoret- 
ical basis for the linear studies in [2], [3], [9] was an abstract approximation framework 
developed by Banks and Ito in [4], It is our purpose here to develop a theory for nonlin- 
ear second order systems that parallels that of [4]. Since we shall use monotone operator 
theory and nonlinear evolution systems (a theory which itself has been shown to possess 
substantial limitations in treating very general nonlinearities), we expect this to be only 
an initial contribution from our studies. However, the theory presented below will allow 
us to treat a large class of nonlinear damping operators that are thought to be of practical 
importance. 

We develop an abstract approximation framework and convergence theory for the identi- 
fication of nonlinear dissipation or damping mechanisms in second order distributed param- 
eter systems. More precisely, we consider the estimation of nonlinear maximal monotone 
damping operators in abstract wave equations with linear symmetric stiffness operators. 
Our treatment here is in the spirit of our earlier work on first order systems in [8] and 
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is an extension of Banks’ and Ito’s results in [4] concerning the identification of second 
order linear systems. By rewriting the underlying second order initial value problem as 
an equivalent first order system in an appropriate product space, we can, in principle, ap- 
ply the abstract approximation theory developed in [8]. However, since the resulting first 
order dynamics are not in general described by a strongly monotone operator, the results 
we present below cannot be obtained via a direct application of our first order theory. 
Only the essential ideas underlying the general abstract Banach space theory involving 
m-accretive operators given in Section 2 of [8] are directly applicable. 


In Section 2 below we review the abstract existence, uniqueness, regularity, and approx- 
imation results for nonlinear evolution equations in Banach space with dynamics governed 
by m-accretive operators. In the third section we consider inverse problems for abstract 
wave equations with nonlinear damping. We reformulate these second order equations as 
equivalent first order systems, and discuss existence and regularity of solutions. The ap- 
proximation theory and convergence results are presented in Section 4. In the concluding 
section we illustrate the application of our theory with an example involving estimation of 
nonlinear velocity dependent damping in a one dimensional wave equation. 


2. Abstract Nonlinear Evolution Equations in Banach Spaces 


Let Xq be a Banach space with norm |-| 0 . Let T > 0 and suppose that Fo G 
Li(0,T]X 0 ). Let A 0 : X 0 — > 2 Xo be an, in general, nonlinear set-valued operator with 
Dom(Ao) = {x G X 0 : f 0}. We assume that for some w G R the operator A 0 + loI 

is m-accretive. That is, that (i) |*i - z 2 |o < K 1 + Aw )( x i ~ x 2 ) + A(yi - y 2 )|o for every 
x u x 2 G Dom(A 0 ),yi G A 0 x u y 2 G A 0 x 2 , and A > 0, and (ii) K(I + \(A 0 + u>I)) = 
U r(E Dom(Ao) (1+ A (Ao +wJ))x = Xo for some A > 0. It is well known (see [10]) that the 
operator A 0 +ujI being m-accretive implies that for each A > 0 the resolvent of A 0 +u>I at 
A, J( A; A 0 +u>I) : X -> X, given by J( A; A 0 +u>J) = (I + A (A) +ul))~ l is a single valued, 
everywhere defined, nonexpansive (nonlinear) operator on Aq- 
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We consider the quasi-autonomous initial value problem 

(2-1) *(<) + A 0 x(t) 3 F 0 (t), 0 < t < T 

(2-2) x(0) = x° 

with x G Dom(.4 0 ). By a (strong) solution to the initial value problem (2.1), (2.2) we 
shall mean a strongly continuous function x : [0, T j -»■ X 0 which is absolutely continuous 
on compact subintervals of (0,T), differentiable almost everywhere on (0,T), and satisfies 
F 0 (<) - x(t) G A 0 x(t) for almost every t G (0,T) and a:(0) = a: 0 . 

A two parameter family of nonlinear operators, { U 0 (t,s ) : 0 < s <t < T}, defined on a 
subset of A 0 is said to be a nonlinear evolution system on i l if it satisfies (i) 

U 0 (t, s)U 0 (s, r)ip = U 0 (t, r )tp for every <p 6 Q and r, s, t with 0<r<s<t<T, and (ii) the 
mapping ( t,s ) — * U 0 (t,s)<p is continuous from the triangle A = {(t, a) :0<cr<r<T}c 
IR 2 into X 0 for each tp G O. 

It can be shown (see [8], [12], and [14]) that with A 0 ,F 0 , and x 0 as they have been 
de fined abov e, there exists a unique nonlinear evolution system {C/ 0 (t,s) : 0 < s < t < T) 
on Dom(Ao) which satisfies 

(i) |[/o(<, •$)<,? - Uo(t,s)ij>\ Q < \<p - ip | 0 , for all <p,ip € Dom(A 0 ) and ( t,s ) G A, 

(ii) \U 0 (s + t,s)<p- U 0 (r + t,r)<p\ 0 <2/ e w(< “ r > |F 0 (r + s) - F 0 (r + r)| dr, for all 

o 

ip G Dom(Ao) and all < > 0 such that s + 1, r + t < T, 

(iii) If the initial value problem (2.1), (2.2) has a strong solution x , then x(t) = Uo(t,s)x(s) 
for all (t,s) G A. 

The strongly continuous function x defined by x(t) = U Q (t,0)x°,t G [0,T], is sometimes 
referred to as the unique mild or generalized solution to (2.1), (2.2). It is immediately 
clear from (iii) above that when the initial value problem (2.1), (2.2) has a strong solution, 
it coincides with the unique mild solution defined in terms of the nonlinear evolution 
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system {U 0 (t,s) : 0 < s < t < T} on Dom(4 0 ). Theorem 2.1 below is the fundamental 
abstract approximation result upon which the approximation and convergence theory for 
second order systems to be developed below will be based. It is similar in spirit to the 
approximation results for nonlinear evolution systems given in [13] and [15]. 

For Banach spaces X and Y we denote the Banach space of continuous linear operators 
from X into Y by C(X, Y). 

Theorem 2.1 For each n G Z + = {1, 2, 3, . . . }, let X n be a Banach space with norm 
| ■ |„. Let tt„ € C (X 0 ,X n ) with | 7 r n y?| n < v \<p\ 0 where v is a positive constant independent 
of n. Let F n G X n ) and let A n : X n -* 2 X " be an operator on Dom(A„) C X n 

with A n +u>I m-accretive. Suppose that there exists a function g € Li(0,T) for which 
|T'„(t)| n < g(t), for almost every t G (0,T) and that 

(i) for each y? € Dom(A 0 ) there exists {^n}^ with <p n € Dom(A„) such that 

lim |y?„ - it n ip\ n - 0, 

n— >oo 

(ii) lim | F n {t) - 7t n F 0 (t)\ n = 0, a.e. t G (0,T), 

n-^oo 

(iii) for some A > 0, we have lim \J (A; A n +u>I)ip n — n n J (A; A 0 +ul)(p\ n = 0 whenever 

v 7 n— »oo 

v>n € X n and <p G Xo with lim \ip n — ^n^>\ n = 0- Then if {U n {t, s) : 0 < s < t < T} is 

n— + od 

the evolution system on Dom(.4„) generated by A n and F n , we have 

(2.3) lim \Un ( t , s) ffn - KnUo ( t , s) (p\ n - 0 

uniformly in s and t for ( t,s ) G A, whenever ip G Dom(Ao) and y? n G Dom(A n ), 

with lim |y? n — 7rn^| n = 0. 

n— ►oo 

The proof of Theorem 2.1, which we omit, is essentially the same as the proof of Theorem 
2.2 in [8]. Some modification, albeit relatively simple and straightforward, is necessary 
however due to the fact that in our application of Theorem 2.1 below we cannot choose the 
In’s as subspaces of Xq. In the proof this is handled by simply expressing all convergence 
conditions and results in terms of the mappings 7r n . (Compare the statements of Theorem 
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2.1 above and Theorem 2.2 in [8]). We note that the convergence result (2.3) is closely 
related to the notion of factor convergence as given in [16]. We remark that this situation 
(-T n (JL -To ) arises frequently in parameter estimation problems and the reader can find a 
more complete discussion with examples along with a linear semigroup version of the above 
approximation theorem in [5, Chapter II]. Whether or not (2.3) implies that U n (t,s)<p n 
and U 0 (t,s)(p are close in any reasonable or useful sense depends of course upon how the 
spaces X n , the norms | • |„, and the mappings 7T„ are chosen. In our treatment below 
we shall be able to apply the following corollary, which follows from Theorem 2.1 via a 
straightforward application of the triangle inequality. 

Corollary 2.1 Suppose that there exists a Banach space X with norm H* for which X 0 
is set-equivalent to X, for each n = 1,2,..., X n is a subspace of X , and all of the norms 
| • |„,n = 0, 1,2, . . . are uniformly equivalent to | • \ x . Suppose further that Iim n n p> = 

n — ► oo 

ip,¥> E X. Then conditions (i) - (iii) in the statement of Theorem 2.1 can be replaced by 

(i) Dom(ylo) C n ^' !X) Dom(d n ) , where Li denotes the lower limit of the sequence of sets, 
{Dom( J 4„)} (see [11], p. 335), 

(ii) Jim °F n (t) = F 0 (t ) in X for a.e. t G (0,T), and 

(iii) for some A > 0, Jim^ J (A; A n -j- ul) <p n = J (A; A n + a >/) <p in X for each <p G X when- 
ever p n G X n with lim <p n = (p in X. Moreover, the conclusion (2.3) can be replaced 
by 

Jim U n (t, s)<p n = U 0 (t,s)(p 

in A , uniformly in t and s for (f , s) G A, for each tp G Dom(^4o) and <p n G Dom(A n ), n = 
1 , 2 ,... with lim p> n = ip in X. 


3. The Identification of Second Order Systems 

Although the identification of nonlinear dissipation mechanisms in infinite dimensional 
second order (mechanical) systems is of primary interest to us here, we treat the somewhat 
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more general problem which includes the simultaneous estimation of linear stiffness. It will 
become clear from our discussions below that an approach similar to the one we use here 
could be used to develop techniques for the estimation of nonlinear stiffness in the presence 
of linear damping. The extension of our general approach to handle inverse problems 
involving the simultaneous identification of both nonlinear stiffness and damping is at 
present still under investigation. We hope to be able to report on results in this direction 
in the not too distant future. 

Let H be a real Hilbert space with inner product < v > and corresponding norm | • ]• 
Let V be a reflexive Banach space with norm denoted by || • ||, and assume that V is 
densely and continuously embedded in H. If we let V* denote the dual of V it follows that 
y o_> }J <_+ V* with the embedding of H in V* dense and continuous. The continuity of 
the embeddings implies the existence of a constant // > 0 for which \tp\ < fi IMI e V * 
and H^ll, < p\<p\,<p G H, where ||.||„ denotes the uniform operator norm on V*. We 
shall also use < •,• > to denote the usual extension of the H inner product to the duality 
pairing between V and V*. Let Q and Z be metric spaces and let Q be a fixed, nonempty, 
compact subset of Q. 

For each q E Q let the operator > 1 ( 9 ) £ C(V,V*) satisfy the following conditions. 

(Al) (Symmetry) For all <£>,?/> £ V, < A(q)<p,ip >—< >■> 

(A2) (Continuity) For each 6 V the map q A(q)<p is continuous from QcQ into V*; 

(A3) (Equi-V-Coercivity) There exist constants «i e R and o 0 > 0, both independent of 
q G Q for which < A(q)<^>,^ > +<*i l^l 2 — a ° llv , l| 2 ) ^ or all yp G F and q G Q\ 

(A4) (Equi-Boundedness) The operators A{q) are uniformly bounded in q for q G Q; that 
is, there exists a constant a > 0, independent of q G Q, for which \\A(q)<p\l < a ||^||, 
for all ip G V. 
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In addition, for each q £ Q let the operator B(q ) : Dom(£(g)) C V — ► 2 V * satisfy the 
following conditions: 

(Bl) (Domain) Dom(6(?)) = Dom(5) is independent of q for q £ Q , and 0 € Dom(5); 

(B2) (Continuity) For each <p £ Dom(5), the map q — * (B(q)<p is lower semi-continuous from 
Q C Q into 2 V (see [17], Section 18.i, page 173); 

(B3) (Maximal Monotonicity) For all (<pi,t/>i),(<p 2 ,fa) € B q = {(^, tf>) £ V x V* : 

P £ Dom (B),^ £ B(q)tp J we have < Vh — ~ ^2 > > 0 with not properly 

contained in any other subset of V x V* for which this monotonicity condition holds; 

(B4) (Equi-Boundedness) The operators B{q) map F-bounded subsets of Dom(5) into subsets 
of V* which are uniformly ^-bounded in q for q £ Q; that is, if S is a F-bounded subset 
of Dom(B), then the set B(q)S is V*-bounded, uniformly in q for q £ Q. 

For each q £ Q let u°(q) £ V, u 1 (q) £ H and / (•; q) £ L x (0, T; H) where T > 0 is 
given and fixed. We shall require that the mappings q —> u°(q ) £ V, q -» u 1 ( q ) g H, and 
q /(<;?), for almost every t £ (0 ,T), are continuous from Q C Q into V,H, and H 
respectively. For each z £ Z, let ;z) be a continuous mapping from (7(0,T; V x II) 
into R + and consider the identification problem given by: 

(ID) Given observations z £ Z , determine parameters q £ Q which minimize the performance 
index 

*%) = <f>((u(q),u(q))-,z) 

where for each q £ Q,u(q) = u(- ;</) is the mild solution to the second order initial 
value problem 

C 3 - 1 ) «(*_)’+ £(«)«(*) + AvM , o < t < t 

(3.2) «(0) = u°(q), u(0) = u 1 ( q ). 
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Following Barbu [10], we rewrite the second order initial value problem (3.1), (3.2) as 
an equivalent first order system in a product space and then apply the abstract theory 
outlined in Section 2 above to make the notion of a mild solution precise. Let X = V x H 
be the Banach space endowed with the norm | • | given by 

1 

KvbVOU = (iMI 2 + l^l 2 ) • 

In addition, for each q G Q let H(q) be the Hilbert space which is set equivalent to X and 
has as inner product < • > q given by 

(3.3) < ,(^ 2 ,^ 2 ) ><?=< ■A(fl , )v ? i>V , 2 > +«i < ¥ , i>¥ , 2 > + < ^iA ’2 > ■ 

Denote the corresponding induced norm by | • \ q - Conditions (Al) - (A4) above guarantee 
that (3.3) indeed defines an inner product on X and that the Banach spaces {' H(q ), M,} 
are uniformly (in q , for q G Q) norm equivalent to the space {X, | • |x}- 

For each q G Q define the operator A(q) :Dom (A(q)) C H{q ) — ► 2 n by 
A (q) = (-Vb {A (5) ¥ + B (q) n H), 

with Dom(A(tf)) = {(</>, V>) G V x V : iff G Dom(B) , {A (q) + B (q) t/>} n H ± 0}- 

Theorem 3.1 There exists an u € R, independent of q G Q, for which the operator 
A(q) +u>I is m-accretive. 

Proof. Our proof is analogous to the one given by Barbu [10] in the case where the 
space and operator do not depend upon a parameter. We first show that for u> G R 
sufficiently large the operator A(q) + u)I is monotone. For (</>,*/>), (f?, 0) € Dom(A(g)) 
and a; G R let (£,$) G (A (q) + ul) (y, V>) and (i ),9) G (A(q) + ujI) (q,0). It follows that 
£ = + utp, rj = + W 77 , = A(q)ip + ip + uip and 6 - A(q)q + 0 + u8 for some 

$ G B(q)rj> and some 9 G B (q) 9. Then, using the monotonicity of the operator B(q) 
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(condition (B 3 )) we obtain 


< (9, VO ~ (rj, #), ( 9 , VO - (v, 9) > g 

=< + 9, A(q) (<p - Tj) + ^ - 0), {tp - T], V> - 9) > q +u; I (<p - T], V» - 0)\ 2 q 

= — ot\ < — 0,ip — rj > — 9 ,tj) — 9 > | (<p — 77 , V’ — 0)| 2 

> -a i < V> “0,9 - »7 > +w < A(q)((p - r)),cp - 7 > +u ;«! |<p - ^| 2 + w j?/> - 6»| 2 

> 0 


if 


a; > wq = sup 


ax < u,v > 


< A (q)u,u > +ai |u| 2 + |u| 2 
We note that u> 0 exists since 


;(«,«)€ X,(u,v) ^ 0, 9 <= Q l . 


a l < 77 , V > < <* 1 // ||u|| |u| 

< A(q)u,v > +01 M 2 + |u | 2 “ a 0 ||u || 2 + |u | 2 

< + M 2 } < aifl 

q-o ||a|| 2 + |u| 2 _ 2min(a 0 ,l)' 

If ^0 < 0, we set u> = 0 and henceforth assume without loss of generality that u> > 0. Since 
A{,q) T tel is monotone, it is also accretive. Indeed, letting A > 0 , we see that 

1 ( 9 , VO - (77, 9 ) \ 2 q < |(v?, 0 ) - (77, 9 ) I 2 

+ A < (9, V>) - (77, ( 9 ), (9, - (77, 0) > ? 

=< (9, VO ~ (*?, 0 ) + A{(9, VO - (>?, #)}, (9, V’) - (77, 0 ) > q 
< 1 ( 9 , V’) - (7 7, *) + A{(& VO - (77, ( 9 )} I, 1(9, VO - (77, 9 ) \ q . 

Thus we obtain 


- Ob 0)U < 1(9, V) - ( 77 , 0) + A{(£, - (rj, 9)}\ g . 

To demonstrate m-accretivity, we let A > 0, and let (77, 9 ) G 'H(q) be fixed but arbitrary, 
and show that (I + A (A (q) + ujI)) (9, xf>) 9 (77,#) for some (9, ip) e Dom^g)). This 
statement can be equivalently written as 

( 3 . 4 ) y? + A (—'ifi + ujip) — rj 
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(3-5) 


ip + X (A(q)<p + B(q)ip + uip) 3 o. 


Solving for <p in (3.4) and the substituting into (3.5), we obtain 

(3.6) </? = ( 1 + Aw) 1 (\‘ip + ri) 

(3.7) (1 + Xu) ip + A 2 (1 + A u)' 1 A (?) ip + XB (?) ip 3 0 - A (1 + Aw) 1 A (?) 7. 

Define the operator T\ (?) G ^(^, V*)by 

(3.8) T x (?) = 0 + Xu)I+ A 2 (1 + Ao;) ! -4(?). 

It follows from condition (A3) that for u 6 V 

(3.9) < T\ (?) u,u > > T\ (A) |u| 2 + t 0 (A) ||ix|| , 

where rj (A) = (1 + Xu) - A 2 (1 + Aw) -1 and r 0 (A) = A 2 (1 + Xu) 1 o 0 . For A > 0 
appropriately chosen (small if ai > 0; any A > 0 if < 0), it is seen that n( A) > 0 and 
hence that T x (?) is monotone and coercive in the sense of Barbu [10, p. 34], That is, 

< T\ (?) u, u > / ||u|| — » 00, as ||u|| -» 00. 

Therefore, applying Barbu’s [10] Corollary 11-13 with his X, B and A chosen as X = 
V,B = T\(q) and A = A £?(?), we find that %{T\{q) + XB(q)) = V*. Consequently there 
exists a ip G Dom (B (?)) for which (3.7) holds. If we then obtain <p from (3.6), it follows 
that H ( I + A (A (?) +ul)) = n (?) and the theorem is proved. 

For each ? G Q define F G Li (0 ,T; W(?)) by F(t;«) = (0, /(<;?)), a.e. t G (0,T), and 
set x° (?) = (u° (?) (?)) G H{q). Theorem 3.1 and the discussions of Section 2 yield 

that A(q) and F(-; ?) generate a nonlinear evolution system {U (<, s; ?) : 0 < s < t < T) on 
Dom(A(?)). Henceforth we shall assume that x°(?) G Dom(A(?)j for each ? G Q, and by a 
mild solution, «(?), to the second order initial value problem (3.1), (3.2) we shall mean the 
V-continuous (recall the uniform norm equivalence of X and 7i(q)) function «(•;?) given by 
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the first component of the 'H(q) (or X) continuous function x ( t ; q) = U {t, 0; q) x° (q) ,t £ 
[0, T 1 ]. We shall take u (•,' ?) to be the ZZ-continuous second component of x(-;?). 


At this point some remarks regarding the closure of the set Dom(A(?)) are in order. We 
show that under a slightly more restrictive boundedness condition on the operators B(q ) 
than condition (B4), the operator A(q) is in fact densely defined; that is, that Dom(A(?)) = 

W( 9 ). 

Theorem 3.2 If condition {BA) is replaced by the stronger condition that the operators 
B{q) map ZZ-bounded subsets of Dom(R) into subsets of V* which are uniformly V*- 
bounded in q for q £ Q, then Dom(A(?)) = ZZ(?). 

Proof: Let (^,0) 6 H (q) and set 

(^n,0„) = j {^\ A{q) + ul) {y,xl>) £ Dorn {A (q)) , 
for n = 1, 2, . . . . As in the proof of Theorem 3.1, we find 
(3-10) <p„ = (l + 1 (n _1 0n + <f) 

C 3 - 11 ) (r n ~! (?) + n -1 R(?)) 0„ (l + n _1 u?) _1 A (q) ip 

where the operator T n - 1 (q) € £ (]/ V*) is given by (3.8) with A = n~ l . Recalling (3.9), 
we have for u £ V 

< T n -i (?) u, u > > T\ (n _1 ) |u| 2 + T 0 (n _1 ) ||u|| 2 , 

where t\ (n ) = (l + n — n 2 (l + n _1 u)) oj > 0 for n sufficiently large, and 

r 0 (n *) = n~ 2 (l + n~ l (x>) 1 q 0 . Note that lim Tj (n _1 ) = 1 and lim r 0 (n _1 ) = 0 

n-*oo v 7 n— »-oo V * 

Now from condition (B3) we find that for some 0 n 6 B (?) 0 n and any 6 £ B (?) (0) 
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(recall that 0 6 Dom (B) by condition (Bl)), 

T i (n _1 ) IV’nl 2 + T 0 ( n_1 ) HV'nH 2 < < ^n- 1 (?) V’mV’n > 

< < T n - 1 (?)</>„, </>„ > +n~ l < 0» > 

=<■!/> — n -1 (l + n -1 u;) * -4 (g) </> — n 1 ^V’n> 

< M IV’nl + ” _1 (1 + ™ -1 ^) 1 ||^(?)¥ ? ll* ll^nll + n 1 11^11*11^ 

< 


) 1 + n 1 
^ wi* + IW 2 + 

+ T -^p- ||0„|| 2 + Ppll^ll 2 + IIWI 2 ■ 


Therefore 

^ + =$p ^ w + i + “gw.- 

It follows that {ip n } is uniformly bounded in H and that {n _1 ||^„||} is also uniformly 
bounded. From (3.11) we obtain 

— 6 —n~ 1 ui^p n — n~ 2 (1 +n-'u)~' A(q)4>n-n~ l (1 + n-'u.)"' A(q) V -n~'B (q) 


Therefore for any xj> n 6 B (q)tp n , we have 

\\ip n - n* < n _1 “> IIV’nIL + n~ 2 (1 + n-'ui) -1 ||./4(?) iML 
+ n -1 (l + n -1 u;) 1 \\A (q) V?||* + n 1 1|^« |1* 

< n _1 w/i \ti> n \ + n~ 2 a ||V>„|| + n~ 1 a ||y>|| + n 1 1| V’rz || * • 

This estimate together with the uniform boundedness of \ip n \ and n 1 \\ip n || and the hypoth- 
esis of the theorem imply that lim ip n = ip m V*. Thus using the uniform boundedness 

n— *oo 

of {\ipn\} and the density of V in H, we obtain ip n -*• weakly in H as n -» oo. 


From (3.10) we find that 


IV’nl < n 1 It/’n I + M 


and 

||V?n|| < HV’nll + M 
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and therefore that the {p n ) are uniformly bounded in both H and V. Since (l + n _1 u>) p n 
= n~ 1 ip n + p, we have 

\pn -<p\< n~ l \ip n \ + n~ l u> \p n \ 

< (n _1 + n~ 2 u) \ip n \ + n~ l bj \p\ — > 0 

as n > oo. Thus lim p n — tp in H. Using the density of H in V* and the uniform 

n—+ oo 

U-boundedness of the tp n , it is not difficult to argue that p n — > <p, weakly in V . Therefore 
{jpn, */’ n ) * {}P, V 7 ) weakly in A and thus the weak closure of Dom is all of X . Since 

A(q) +ul is m-accretive by Theorem 3.1, the strong closure of Dom { A{q )) is convex. It 
follows, therefore, that Dom(j4(5)) = 'H(q) and the theorem is proved. 

When we consider examples in Section 5 below, we shall show that there is an im- 
portant class of nonlinear models for dissipation which satisfy the boundedness condition 
hypothesized in the statement of Theorem 3.2. 

Finally, we note that another class of operators B(q) which lead to the operator A{q) 
being densely defined is the one which includes operators B{q ) 6 £{V,V*) for which the 
set {p 6 V : B(q)<p 6 H } is dense in H. For then condition (A3) implies that {p € V : 
H } is dense in V (see [22], Theorem III.2.B) and consequently Dom(A(^)) = 
{(¥b ip) £V x V : A(q)p + B (q) ip € H } is dense in 7i(q). In particular, this will be the 
case for the well known linear Kelvin-Voigt (i.e., B{q) ~ A (q)) and the so-called structural 
(i.e. B(q) ~ A(q)? ) viscoelastic damping models. 

4. Approximation Theory and Convergence Results 

We apply Ritz-Galerkin finite element techniques to discretize the state equation. For 
each n = 1,2,... let H n be a finite dimensional subspace of H with H n C V. Let 
P n : H — » H n denote the orthogonal projection of H onto H n computed with respect to 
the H inner product. We assume that jP n Dom(£?) C Dom(Z?) and that the approximation 
condition 

(P) lim || P n p - p\\ = 0, p G V, 

n — >oo 
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is satisfied. We note that condition ( P ) implies that lim \P n ip - <p\ = 0 , tp G H as well. 

For each n = 1,2, ... and q £ Q, let A n (q) € C(H n ) = C(H n ,H n ) and B n (q) .Dom 
(5„) C H n 2 Hn denote the usual Galerkin approximations to the operators .4(g) and 
B(q) respectively. More precisely, for ip n 6 H n let A n (q)q>n = V’n where is the unique 
element in H n guaranteed to exist by the Riesz representation theorem satisfying 
< A(q)<Pn,Xn > = < >,Xn £ H n . Similarly, for <p n EDom (B n ) = Dom(B)nH n ? 

0, let B n (q)ifin = {V> '» € H n :< </>,X« >=< *Pn,Xn >,X» € H n , for some rf> e B (q) <p n ] • 
We set f n (-,q) = P n f(",q ) 6 Xj (O.TjJT.), «£(«) = w n(<0 = 

Our approximation framework is based upon the following sequence of approximating 
identification problems. 

(ID„) For the observations z £ Z given in problem (ID), determine parameters q n E Q which 

minimize the functional 

$»(?) = <P( u n(q),Un(qy,z) 

where for each q £ Q u n (q ) = «„(•; ?) is the mild solution to the second order initial 

value problem in H n given by 

( 4 . 1 ) u n (t) + B n (q)un(t) + An(q)un(t)Bfn(t-,q)i 0 <t<T, 

(4.2) u n (0) = u° n (q), «„(0) = v} n (q). 


Once again to establish the existence of a unique mild solution for each q E Q, and to 
develop a convergence theory, we reformulate the initial value problem (4.1), (4.2) as an 
equivalent first order vector system in an appropriate product space, and then apply the 
abstract existence and approximation results outlined in Section 2. 

For each n = 1,2,... let X n = H n x H n and for each q eQ let H n (q) denote the space 
X n considered as a subspace of the Hilbert space H(q). Define the operator A n (q ) : Dom 
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(4 n (g)) C H n (q) -► 2 W »<«> by 

A-n(q) (V^n? V’n) ( ^ni A n (q)(p n ~b ^n(9)V , n) 

for each (^ n ,ipn) € Dom(,4 n (g)) = H n x Dom(5 n ). Set F„(f;g) = (0,/„(t;g)) for almost 
every t G (0, T) and let 4(g) = (u° n (g), 4(g)). Note that F„(-;g ) G Ti (0, T; ft„(g)) and 
that the assumptions that z°(g) G Dom(.4(g)) and P n Dom(5) C Dom ( B ) imply that 
4(«) ^ Dom(d( ? )). 

With the above definitions of ,4„(g) and #„(g), it is not difficult to argue that A n (q)-\-uI 
is m-accretive in H n (q). (The proof is analogous to the proof of Theorem 3.1.) It follows 
that for each n = 1, 2, . . . and each q G Q, A n (q) and F n (-] q ) generate a nonlinear evolution 
system, {U n (t, s; q) : 0 < s < t < T} on Dom(.4(g)). By a mild solution u n (g) — u n (-; g) 
to the initial value problem (4.1), (4.2) we shall mean the F-continuous first component of 
x n{-',q ) = U n (-,0; g) x°(g). The if-continuous second component of x„(-;g) will be taken 
to be u n (•; g). 

The primary result of this paper is given in the following theorem. 

Theorem 4.1 . For each n = 1,2,... the problem ( ID n ) admits a solution q n G Q. 
Moreover, the sequence {gnJ^lj has a convergent subsequence {g n *}^Li with lim q Uk = 

fc—>oo 

q £ Q) where q is a solution to problem (ID). 

Proof: In light of the compactness of Q and the continuity of <f>, the existence of a solution 

q n £ Q to problem ( ID n ) can be established by demonstrating that the mapping q — ► 
x n( m ; q) is continuous from Q C Q into C(0, T; X ). Furthermore, as we have demonstrated 
in several of our earlier papers (see, for example, [6], [8]), showing that q — lim q nk is 

k—+oo 

indeed a solution to problem ( ID ) requires arguing that for any convergent sequence 
with Qn G Q and lim g„ = g 0 G Q we have 

77- — ► OO 

(4.3) Jin ^ \ x n (^j qn) x (^j 90 )^ = 0 


16 


uniformly in t for t £ [0 ,T]. Both the continuous dependence and the convergence (4.3) 
can be demonstrated via an application of Corollary 2.1. We shall establish (4.3) here; 
continuous dependence can be argued in an analogous manner. 

Let C Q with lim q n = q 0 in Q. The compactness of Q of course implies that 

q 0 £ Q. In applying Theorem 2.1, we take X 0 = 'H(q 0 ),X n = ^n(?n), H„ = M ?n = 
0 , 1 , 2 ,... , and define 7 r„ £ C{'H{qo) ,'H n {qn)) by -Kn^^xp) = (PnVh Pn*f>) , (v 7 ; V > ) £ 

(go)- Note that condition (P) implies that the 7 r n are uniformly bounded and converge 
strongly in X to the identity. We set F n = F n q n ) , Fq = F(- m ,qo),A n = A n (q n ), Aq = 
A(q 0 ), and will use Corollary 2.1 to establish that 

lim | U n (t,0-,q n )x° n ( q n ) - U (t,0-,q o )x° (?o)| x = 0 

71—^00 1 

uniformly in f, for t £ [ 0 , T]. Under the general assumptions that we have made above, all 
of the hypotheses of Corollary 2.1 can be verified immediately with the exception of the 
resolvent convergence condition (ii). We establish it here. 

Let A > 0 be fixed, let (<p,ip) £ 'H(qo), (v’n.V’n) £ H n (q n ), and assume that 
lim \((p n ,ip n ) — (v^VOIx = 0- It then follows that lim = ip in V and lim^ip„ = ip in 
H. Setting ( 77 , #) = J(A; A(qo) + ul)((p,ip) £ Dom (A(g 0 )) an d 
(r?„, 9 n ) = J (A; A n (q n ) + ul) (v?„, ip n ) € Dom(A„ (?„)) we obtain 

( 4 . 4 ) q = (1 + Aw) ^Afl + ^j) 

( 4 . 5 ) (T\ (q 0 ) + \B(qo))0 3xf>-\(l + A u>)~ 1 A (go) 

(4.6) T/n — (1 + Au;) (A 9 n + (p n ) 

(4.7) (?n) + AZ? n (q n )) 3 n A (1 + Au>) ./4n (?n) 

where T\(q) is given by (3.8), and 

T nt \ ( q ) = (1 + Au>) / T A 2 (1 + Ao;) A. n (?) £ F(H n ). 


17 


Recalling (3.9), we choose A > 0 so that tj(A) > 0. We first show that the 8 n £ H n 
are uniformly V-bounded in n. From (3.9), conditions (A4), (B3), and (B4), and the 
definitions of the operators A n (q) and B n (q ) we find for some 9 n 6 B n (q n )0 n and any 
in £ B(q n ) (0) 


r 1 (A)|^| 2 + r 0 (A)||^|| 2 < < T x (q n )8 n ,9 n > 
= < T ni x(q n )0 n ,6 n > 

- < < Tn t \(q n )9n,9n > +A < 8 n -i„,0 n > 


— < ‘ipn — A(l + Au>) 1 A(q n )<Pn — A£„,#„ > 

< |0»| \0 n \ + A(1 + Aw)" 1 a \\<p n \\ \\9 n \\ + A/? 0 \\9 n \\ 

< -4- |V-„I 2 + 4^ K\ 2 + “ 2(1 + Ai4,r ‘ 


+ 


2t 1 (A) 
To (A) 


<*0 


II*. 


2 ,( 1 + Aw) 2 T 0 (A) 


+ 


«0 




IIM S 


where /? 0 = sup |||£||* : i G B (?)( 0),9 6 Q^. It follows that 


T,(A) i«„i 2 + ^im 2 < 


1 


IW J + - (1 + Aw) "‘ 


IMI 2 + (1 + 


2 2 -2r,(A) ' «„ ” «„ 

The fact that lim \\(p n — q>\\ — 0 and lim \jp n — ip\ = 0 implies that ||cz> n || and \rb n \ are 

n — *o© n — *oo 11 ' 1 

uniformly bounded, and consequently that the 0 n are uniformly V-bounded, 


Once again from (3.9), for any 6 n (q n ) G B ( q n )6 n , 9{q n ) G B(q n )9, and d(q 0 ) G B(q Q )6 
we have 


n (A) \9n - 9 \ 2 + To (A) || 9 n - 9 \\ 2 < < T A ( 9n ) {<? n - 9} ,8 n - 0 > 

< < T x {q n ) {#n - 6} A - 0 > +A < 9 n (q n ) - 9(q n ),0n ~0> 

= < T n>A ( 9n ) 0 n - T A ( 9o ) > +A < 9 n (q n ) -0(q o ),0 n - p n 8 > 

+ <T X (q n ) 9 n - T x ( 9o ) 0, P n 9-9>+\< 8 n (q n ) -8{qo) ,P n 0-0> 

+ < T x (q o )0 - T x (q n )9,0 n - 8 > + A < 8(q 0 ) - 8{q n ),0n — 9>. 


Now from (4.7) there exists a 0„ € B n (q n )9 n for which 


Tn,\ ( 9n ) + A#„ — — A (1 + Aw) 1 A n ( 9 n) 
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Choose 9 n (q n ) to be any element in B(q n )9 n for which 


^ 8n (?n) > Xn ^ ^ Xn > Xn ^ f^n - 

That such a ^ n (g„) € B(q n )8 n exists is “guaranteed” by the definition of B n (q n )- Also, 
recalling (4.5), we choose 9 (qo) E B(qo)8 so that Ta(<Jo)0 + A 9(qo ) = 

%l> — A (1 + Au>) 1 A (qo) <p- Then 

t\ (A) \6 n — 8\ 2 + r 0 (A) \\9 n -9\\ 2 

< < r n , A (q n )6 n + \8 n - {T a ( 9o ) 9 + A 9 (q 0 )},8 n - P n 8 > 

+ < T a (5„) 9 n - T x (q 0 ) 8, P n 9-9>+ A < 8 n (q n ) - 9 (q 0 ) ,P n 9 — 8 > 

+ <T x {q 0 )8-T x (q n )9,9 n -9 > +A < 9(q 0 ) - 9(q n ) ,9 n - 9 > 

= < i>n — 1p,8 n - P n 8 > +A (1 + Aa;) -1 < . 4 ( 50 ) V 3 - A n (<?n) — P n 9 > 

+ <T X (q n ) 8 n - T x ( 50 ) 8 , P n 8-8>+ A < 9 n (q„) - 9 ( 50 ) ,P n 8-9> 

+ <T X (q 0 )8- T x (q n )9,9 n -9>+ A < 8(q 0 ) - 8 (q n ) ,8 n -9> 

= < tf>n — V’j 8 n — P n 9 > +A (1 + Au;) < A ( 50 ) y ~ A ( q n ) V 3 ? 8 n ~ P n8 > 

+ A (1 + Au;) -1 < *4 (g n ) (y? n — <p) ,9 n — P n 9 > + (1 + Aw) < 9 n - 8, P n 9 - 9 > 

+ A 2 (1 + Au;) -1 < A(q n ) 9 n -A(q 0 )9, P n 8 - 9 > + A < 9 n (q n ) - 6(q 0 ) ,P n 9-9> 

+ A 2 (1 + Au ;) -1 < (4 ( 50 ) 8- A (q n ) 9, 9 n - 8 > +A < 9 (q 0 ) - 0 (q n ) ,9 n -9> 

< |0 W - V’l |0» - Pn9\ + A(l + Au;) -1 \\A(q 0 )<p ~ A(q n )<p\l || 9 n - PJ\\ 

+ A (1 + Au;) ||«4 ( q n ) i}Pn ~ V 3 )!!* Il^n — -Pn$|| T (1 + Au>) 1 9 n — $| \P n 9 — 9\ 

+ A 2 (1 + Au,)- 1 || A ( ? „)#„ - A (<ro)«||. W - #11 + A||«» (?„) - #(«o)||,||P»# - #11 
+ A 2 (1 + Au.)- 1 \\A (q„) # - A (,„) #||. ||«. - #|| + A||«(*,) -#(}„) ||. ||». - #11 
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5: {|#n| + |^|} IV’n ~ + A (1 + Aw) 1 {||$n|| + K ||^||} ||«4(9o) V ~ A (<Zn) <p||* 

+ A (1 + Aw) _1 a {||0„|| + K ||0||} II Vn - ¥>|| + (1 + Aw ) {|0 n | + 1^1} | P n9 — Q\ 

+ A 2 (1 + Aw)" 1 a {\\6 n \\ + ||0||} \\P n 0 - e\\ + A{||0„ (?„)||. + ||0(<7o)||.} \\PnO - 0|| 

+ A 2 (1 + Aw)- 1 {||* n || + \m \\A(qo) o -A(q n )0\l + A {p„|i + \\0\\} ||*(*>) - 0(?n)||. 

where K is the uniform bound on the operators P n G C(V, V) guaranteed to exist by 
condition (P). By condition (B2) we can choose 9(q n ) G t3(q n )9 so that J.im^0(q„) = 9 (q 0 ) 
in V*. Thus, the fact that lim xf> n = %J) in H , lim <p n = ip in V, the uniform boundedness 

n— ► oo »— *■ oo 

°f ||0 n || (and therefore |0„| as well), the continuity condition (A2), and the boundedness 
condition (B4) imply that lim 9 n — 9 in both H and V . From (4.4) and (4.6) it then 

n— ►oo 

immediately follows that lim rj n = rj in V and the theorem is proved. 

n— * oo 

Remark In practice it is frequently the case that the admissible parameter set Q is 
also infinite dimensional; for example, when the unknown parameters to be identified are 
elements in a function space. In this situation the set Q must also be discretized and 
Theorem 4.1 must, to a certain degree, be modified. For each m — 1,2, . . . let I™ : Q C 
Q _> Q be a continuous map with range Q™ = I™(Q ) in a finite dimensional subspace 
of Q and with the property that lim I™(q ) = q , uniformly in q for q G Q. We then 

771— +00 

define the doubly indexed sequence of approximating identification problems (ID™) by 
letting problem (ID™) be the problem (ID n ) with the set Q replaced by the set Q™ . The 
modification to Theorem 4.1 would state that each of these problems admits a solution 
q™ G Q™ and that the sequence (g” 1 } will have a Q-convergent subsequence whose limit 
is in Q and is a solution to problem (ID). We note that each of the problems (ID™) 
involves the minimization of a continuous functional over a compact subset of Euclidean 
space subject to finite dimensional state space (ODE) constraints. As such, for each m and 
n, problem (ID™) may be solved using standard computational algorithms and techniques. 


20 



5. An Example 


In order to illustrate the application of the general theory that was presented above, we 
consider the problem of estimating or identifying the nonlinear damping term in the forced 
one dimensional wave equation 

( 5 . 1 ) + o<,<i 

with the Dirichlet boundary conditions 

(5.2) tx(t,0) = 0, u(t,l) = O s t> 0 

and initial data 

du 

( 5 . 3 ) u ( 0 , 77 ) = u°(rj), — ( 0 , r?) = u\r]), 0 < r? < 1. 

We note that for definiteness we have chosen the Dirichlet boundary conditions (5.2), how- 
ever all of the discussion below would remain with any of the usual self-adjoint boundary 
conditions (i.e. Neumann, Robin, etc.). We assume that a > 0,/ € Z< 2 (( 0 , T) !x (0, 1)), u° £ 
Hq( 0, 1) and u 1 £ L 2 (0 , 1). Let Q be a metric space and let Q C Q be compact. For each 
q £ Q we assume that the mapping &(§)(•) : R — + 2 R satisfies the following conditions: 

(bl) The set {( 7 ? 6 Hq (0, 1): there exists ai/> £ H~ l (0, 1) such that tp (q) 6 b(q ) (v?(^))> a - e - 
rj G ( 0 , 1 )} is independent of q € Q, and 0 £ 6 (g)( 0 ), 

(b2) The mapping q —y 6 (g)(^) is lower semi-continuous from Q C Q into 2 R for almost every 

6eR, 

(b 3 ) For each q £ Q the mapping &(?)(•) is nondecreasing and for some A > 0 the inclusion 
6 + A 6 (g) (8) 3 f has a solution 8 £ R for each £ £ R (i.e. 6 (g)(-) is maximal monotone 
in R), 

(b 4 ) There exists a polynomial p, independent of q £ Q for which |f?| < p(|0|) for all 8 £ 
6 (g) (| 0 |), and almost every 8 £ IR. 
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To reformulate the initial-boundary value problem (5.1) - (5.3) as an abstract system of 
the form (3.1), ( 3 . 2 ), we set H ~ L 2 ( 0 , 1 ) endowed with the standard inner product, and 

l 

let V = Hq ( 0 , 1 ) be endowed with the norm \\<p\\ = f \D<p(r})\ dr}. In this case we have 

o 

V * = H~ 1 (0, 1 ) together with the dense and continuous embeddings V H <-* V*. We 
define the operator A 6 £(V,V*) by 

l 

(5.4) < >= a J D<p (t)) Dip (t}) dr), p, $ e H] (0, 1 ) • 

0 

It is easily shown and well known that the operator A given by (5.4) satisfies conditions 
(Al) - (A4). (Note that in the present example the operator A is assumed known and as 
such does not depend on g.) 

To define the operator B{q) we follow the treatment in [18] and use the notion of a 
subdifferential of a proper convex lower semicontinuous mapping (see [10]). For each q 
G Q let & 0 (< 7 )(‘) denote the minimal section of the mapping 6 (g)(*). That is &o( 9 )(*) is the 
single-valued mapping from R into R defined by 6 o(g ) = 9 , where 9 is the unique element 
in b(q)(9) of minimal absolute value. That the minimal section of 6 (g)( *) is well defined is 
a consequence of condition (b3). Since Dom( 6 (g)(*)) = R, the proper, convex, lower semi 
-continuous function j(q q) : R — > R can be defined by 

9 

;(#;«)= /«*(«)({)<*£ 

0 

where R denotes the extended real numbers. It can be shown that j(q q) is bounded below 
by an affine function so that j(9;q) > — oo, 9 £ R, and that dj (qg) = 6 (g) (•), where d 
denotes the subdifferential operator. For each q € Q define 7 (•; g) : Hi (0, 1 ) — ► R by 

( , x / j(<p (■);?) e L i (0,1) 

7 (Vb Q) — \ 0 

( +00 otherwise. 

Then 7 (-; q) is also proper, convex and lower semi-continuous and we define B(q) : 

Dom({?) C V — * 2^ by 

(5.5) B(q)v> = dy(v,q). 
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ill Hill' ' 



where p G Dom (B(q)) = {p G H] (0, 1) : d^q) ± 0}. (The fact that Dom (0(g)) = 
Dom(£) is independent of q is a consequence of condition (bl). ) It can be shown (see, for 
example, [10, p. 61]) that for tp G Dom (£),</> G B(q)p if and only if *l>(q) e Ktiiv (*))» 
a.e. q G (0, 1), and via conditions (bl) -(b4), that the operators B(q) given by (5.5) satisfy 

conditions (Bl) - (B4). 

To illustrate the formulation of an inverse problem we take the observation space Z to 
be given by Z = x {R f x I 2 (0, 1)} and the performance index from C(0,T; V x H ) into 
the nonnegative reals to be given by the weighted least-squares functional 

j |u (t t ,q) - z* (q)\* dq) 
i=l j= 1 0 

for (u,u) G C{ 0 ,T‘,V x H) and z = ((*}, z\) , • • • (4, 4)) G ^ with Pi ,a t > 0,z = 
1,2 ,... ,v , 0 < ti < t 2 < ... < t v < T, and 0 < qi < q 2 < m < 1 - 

We provide two specific examples of possible choices of the parameter space Q, the 
admissible parameter set Q, and mappings %)(•) : R -» 2" satifymg conditions (bl) - 
(b4). A relatively simple example might involve the estimation of the constant parameters 
q = (a 0 ,/? 0 ,6>o) in the saturating polynomial function 

_ / A> \0\" s 9 n ( 6 ) » -0o<0<9 o 

^ q ^ -\t3 0 \6 0 \ O0 sgn(6), |0| > ^o- 


In this case we would have Q — R 3 and Q — {(oojAh^o) '• 0 < no < d, 0 < /?o — Pi® — 
0 o < 6} for some a, P, 6 > 0 given and fixed. We note that when o 0 = 0, in order to satisfy 

the maximality in condition (b3) we must take 


Kim = 


/3 0 sgn(6), 

[—3o,Po] 


e^o 

6 = 0 . 


We note also that when a < 2 we have Dom(B) = V and that the operators B(q) given 
by (5.5) map ^-bounded subsets of 7 into ^-bounded subsets, uniformly in q for q G Q. 
Thus the hypotheses of Theorem 3.2 are satisfied and we see that Dom(A(^)) = V X H . 
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A second, somewhat more challenging problem might involve the estimation of a func- 
tional parameter. In this case we take Q — C(g(R), the set of all bounded, continuous, real 
valued functions on IR endowed with the supremum norm. We take Q to be the Q-closure 
of the set 

{qeQ :q(0) = -q(-0) ,0q(0) > 0,for OgR, 

(5-6) | q («)| =q(0o), for \0\ > 0 O , q E H 1 (~0 O ,0 O ) , 

- K °’ v' (&) - f° r a - e ■ 0 € (-00, 0o)} 

where 6 0 and K 0 are given positive constants. It is not difficult to argue (via the Arzela- 
Ascoli Theorem) that Q is a compact subset of Q. For q E Q we set b(q)(6) = q(8),0 6 R. 


With regard to illustrating the approximation or discretization of the state, we briefly 
outline a linear spline-based scheme. For each n = 1, 2, . . . , let H n = span {‘p } n } n ~* where 
( l ) h is th e j — th linear 5-spline on the interval [0, 1] defined with respect to the uniform 
mesh {0, l/n,2/n, . . . , 1}. That is, 


<Pn Of) = < 


( 0 0<J7<i^i 

nrj-j + 1 < rj < { 

j + 1 - nr, i < n < i±l 

- 0 ^ < V < 1, 

j = 1,2,... ,n - 1. It is immediately clear that H n C V = H£( 0,1 ), for each n. Let 

P n . H ► H n denote the orthogonal projection of 52(0, 1) onto H n computed with respect 

to the usual L 2 inner product. Using well known estimates for interpolatory splines (for 

example, those found in [23]) it can be argued that condition (P) is satisfied. 


In the example involving the identification of a functional parameter, the admissible 
parameter set Q chosen to be the (2-closure of the set given in (5.6) can be discretized as 
follows. For each m = 1,2,... and 0 € R+ let (-,0)}™ 0 denote the standard linear 
5-splines on the interval [0,0] defined with respect to the uniform mesh (0 — 2© 0) 

and then extended to a continuous function on the entire positive real line via ipj 1 0^ = 

V>r(0;0M>0. For q (z Q set 

m 

( 5 - 7 ) ( 0) = ^q (^) (\9\\0 q )sgn(0) 

3 = 1 
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for 6 £ R where 6 q is the number in (0, 6 0 ) for which | q (0)| = q (Q q ) , |#| > 0 q . (Note that 
the lower limit in the sum in (5.7) is 1 rather than 0 since q £ Q implies g(0) = 0.) Using 
the Peano Kernel Theorem (see [21] p. 22) it can be argued that 

\\i m t - ?IL = su p ic« - «l < I (£)' /2 K ° 

R L 

and consequently that lim I m q = 5 , uniformly in q for q E Q- We would then set 

m — *oo 

Q™ = I m (Q). 

For H n and Q™ as defined above, the finite dimensional initial value problem (4.1), (4.2) 
takes the form 

(5.8) M n w n (t ) + C n (w n (t)\ q m ) + A n te„(t) = F n (t ), 0 < t < T, 


(5.9) 


M n w n ( 0) = w° n , M n w n { 0) = icj, 


where F n (t),u>° , and are the (n-l)-vectors whose i-th components are given by F\(t) 
= < f(t),<p’ n >,w°' =< u°,(p’ n > and >, respectively, and M n and K n 

are the (n - 1) x (n - 1) matrices whose (ij)-th entries are given by M' n ' j =< <f J n > 
and Rtf =< aq>’ n ,<p J n >. The vector function C n ( ■ ;q m ) : R n_1 -» R n_1 is given by 


Tl 

C' n (v,q™) = J {nx — * + 1} q m ({nx — i) {u* — u ' -1 } + v l * ) dx 


l-l 

n 


+ 


7T 

J {i + 1 — nx} q m ({nx — i} {v‘ +1 — u'} + u*) dx, 


i = 1,2, — ,n - 1, for v £ R n_1 with u°,u n = 0. If u>„(- \q m ) is the solution to 
the second order initial value problem (5.8), (5.9) corresponding to q m £ Q™ , then 
Un(t’,q m )= Y w J n (t;q m )<p J n and u n (t-,q m ) = Y, w 3 n (t;q m )<p{, for t £ [0,T]. If q m £ Q™ 


i=i 


is given by q™(6) = Y Q?'l>Y(\8\\0q tn ) s 9 n ( 9 )> 0 e the identification problem (ID™) 

3 = 1 
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becomes one of determining parameters (jJ", q™,9 q ™) in some compact subset of R m+1 
which minimize $ n (q m ) = <f>{u n {q m ), u n (q m ); z). 

A discussion of some implementation questions relevant to the schemes outlined (in 
particular with regard to supercomputing) together with numerical results can be found 
in [6]. 
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